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Abstract : We advance the novel hypothesis that water in the presence of finite 
localized strain within fault gouges may lead to the phase transformation of stable 
minerals into metastable polymorphs of higher free energy density. Under increas- 
ing strain, the transformed minerals eventually become unstable, as shown from an 
application of Landau theory of structural phase transitions. We propose that this 
instability leads to an explosive transformation, creating a slightly supersonic shock 
wave propagating along the altered fault core leaving a wake of shaking fragments. 
As long as the resulting high-frequency acoustic waves remain of sufficient amplitude 
to lead to a fluidization of the fault core, the fault is unlocked and free to slip under 
the effect of the tectonic stress, thus releasing the elastic part of the stored energy. 
We briefly discuss observations that could be understood within this theory. 



in EARTHQUAKE THERMODYNAMICS AND PHASE TRANSFORMA- 
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1 Introduction 



An earthquake is a sudden rupture in the earth's crust or mantle caused by tectonic 
stress. This premise is usually elaborated by models that attempt to account for 
seismological and geological data as well as constraints from laboratory experiments. 
Notwithstanding almost a century of research since the standard rebound theory of 
earthquakes was formulated [Reid, 1910], the complex nature and many facets of 
earthquake phenomenology still escape our full understanding and is reflected in the 
disappointing progresses on earthquake prediction [Keilis-Borok, 1990; Mogi, 1995; 
Kanamori, 1996; Geller et al, 1997]. 

We first review three important paradoxes, the strain paradox, the stress paradox 
and the heat flow paradox, that are difficult to account for in the present stage of un- 
derstanding of the earthquake processes, either individually or when taken together. 
Resolutions of these paradoxes usually call for additional assumptions on the nature 
of the rupture process (such as novel modes of deformations and ruptures) prior to 
and/or during an earthquake, on the nature of the fault and on the effect of trapped 
fluids within the crust at seismogenic depths. 

We then review the evidence for the essential importance of water and its interac- 
tion with the modes of deformations [Kirby, 1984; Hickman et al, 1995; Evans and 
Chester, 1995; Thurber et al, 1997]. We see then present our scenario. 

2 Strain, stress and heat flow paradoxes 
2.1 The strain paradox 

Within the elastic rebound theory, the order of magnitude of the horizontal width 
over which a shear strain develops progressively across a fault prior to an earthquake 
is of the order the thickness 15 km of the seismogenic zone. This estimate is robust 
with respect to the many ways with which one can refine the model [Turcotte and 
Schubert, 1982; Sornette, 1998b]. 

Modern geodetic measurements are sufficiently precise to test for the existence of 
strain localization. The surprise is that there is no geodetic evidence of the existence 
of a concentration of shear deformation at the scale of 1 — 30 km around major active 
faults in many situations, even if others exhibit it [Pearson et al, 1995]. The geodetic 
displacement profiles obtained across section of faults that have not ruptured in the 
last decades give an essentially uniform strain over distances of 150km or more from 
the fault [Walcott et al, 1978; 1979; Shen et al, 1996; Snay et al, 1996]. 

On the other hand, strain concentrations are observed very neatly after a large 
earthquake [Shen et al, 1994; Massonnet et al, 1994; 1996]. The observed postseis- 
mic relaxation, with concentration in the vicinity of the fault, is attributed either to 
viscoelastic relaxation in the volume of the crust and upper mantle or to afterslip or 
continued slip on the fault rupture planes. 
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2.2 The stress paradox 



2.2.1 Statement of the paradox 

There is a large body of literature documenting the maximum shear stress necessary 
to initiate sliding as a function of normal stress for a variety of rock types. The best 
linear fit defines a maximum coefficient of static friction f s = 0.85 [Byerlee, 1977], 
with a range maybe between 0.6 to 1 approximately []. The range has been confirmed 
by stress inversion of small faults slip data in situ with result 0.6 ±0.4 [Sibson, 1994]. 

There is a big discrepancy between the stress threshold of 200 MPa at 10 km 
depth implied by this value of friction measured for rocks in the laboratory and the 
stress that is available in nature to trigger an earthquake. To develop such a large 
stress, shear strain should reach e ~ 3 10~ 3 [Sornette, 1998b] demanding an intense 
localization over an horizontal width of a few kilometers. As already discussed, this 
is not usually observed. 

In addition, the stress of the San Andreas fault zone deduced from a variety of 
techniques is found to be low and close to perpendicular to the fault [Zoback et al, 
1987]. This is in contradiction with the frictional sliding model. A more detailed 
tensorial analysis assuming hydrostatic pore pressure gives the depth-averaged shear 
strength of faults in the brittle continental crust under a typical continental geotherm 
and strain rate of about 35 MPa in normal faulting, 150 MPa in thrust faulting and 
60 MPa in strike-slip faulting [Hickman, 1991]. A set of stress-field indicators, in- 
cluding borehole breakouts, earthquake focal mechanisms and hydraulic fractures, 
suggests that many (but not all) active faults are sliding in response to very low 
levels of shear stress [Zoback, 1992a; 1992b]. Also noteworthy is the paradox of large 
overthrusts [Brune et al., 1993], namely that thrust faults exhibit an orientation too 
close to the horizontal to obey the usual friction law, thus also requiring an anoma- 
lously small resistance to friction. Such low-angle faulting is observed in many places 
and is in contradiction with Anderson's classification of faults [Anderson, 1951] and 
with the friction theory usually used for the other modes of faulting. The conclusion 
on this weakness of faults is contrary to traditional views of fault strength based upon 
laboratory experiments and creates a serious problem as one cannot rely directly on 
the knowledge accumulated in the laboratory. 

The situation is made even more confusing when one refers to the occasional ob- 
servations of very high stress drops (30 MPa to 200 MPa) in moderate earthquakes 
[Kanamori, 1994]. This would indicate that the stress drop (and therefore the ab- 
solute stress level) can be very large over a scale of a few kilometers. But then we 
should see large strain and anomalous heat flux due to frictional heating (see below). 



1 Heuze [1983] reviews a large body of the literature until that time on the various properties of 
granitic rocks. 
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2.2.2 Proposed resolutions 

To account for these puzzles, many suggestions have been made. Let us mention 
in particular that Chester [1995] has proposed a multi-mechanism friction model 
including cataclitic flow, localized slip and solution transfer assisted friction in order 
to describe the mechanical behavior of the transitional regime at midcrustal depths. 
Blanpied et al. [1995] also use the multi-mechanism friction model for frictional 
slip of granite. They stress that extrapolating the laboratory results to conditions 
not encompassed by the data set (i.e. to approach the conditions in the crust) is 
uncontrolled as many mechanisms are competing in a complex way. The solution 
which is often proposed to resolve this problem is to invoke fluid pore pressure. The 
presence of a fluid decreases the normal stress and thus the shear stress necessary to 
reach the threshold while the most favorable slipping direction is unchanged. This last 
condition ensures the compatibility with the in-situ stress inversion measurements. 
The problem however remains to find a mechanism for pressurizing fault fluids above 
the hydrostatic value towards the lithostatic value in short time scales compatible 
with earthquake cycles [Lachenbruch, 1980; Rice, 1992; Blanpied et al, 1992; 1995; 
Scholz, 1990; Brace, 1980; Lockner and Byerlee, 1993]. To summarize, it seems that 
invoking fluid overpressure to weaken the effective friction coefficient serves a single 
purpose and creates novel difficulties to interpret other observations. We refer to 
[Sornette, 1998b] and the references therein for more information. 

2.3 The heat flux paradox 

2.3.1 Statement of the paradox 

The heat flow paradox [Hickman, 1991] in a seismically active region was first pro- 
posed by Bullard [1954] : to allow for large earthquakes, a fault should have a large 
friction coefficient so that it can store a large amount of elastic energy and overpass 
large barriers. However, if the dynamical friction coefficient is large, large earthquakes 
should generate a large quantity of heat not easily dissipated in a relatively insulating 
earth. Under repetition of earthquakes, the heat should accumulate and either result 
in localized melting (which should inhibit the occurrence of further earthquakes) or 
develop a high heat flow at the surface. Observations over the entire state of California 
have shown the absence of anomalous heat flow across the major faults [Henyey and 
Wasserburg, 1971; Lachenbruch and Sass, 1980; 1988; 1992; 1995; Sass et al, 1992]. 

2.3.2 Proposed resolutions 

A standard explanation is that no significant heat is generated because the dynam- 
ical friction coefficient is low on the most rapidly slipping faults in California. Two 
classes of models are then usually proposed. In the first class, the low friction is 
produced dynamically during the event itself. Various mechanisms are invoked, such 
as crack-opening modes of slip [Brune et al, 1993; Anooshehpoor and Brune, 1994; 
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Schallamach, 1971], dynamical collision effects [Lomnitz-Adler, 1991; Mora and Place, 
1994; Pisarenko and Mora, 1994; Schmittbuhl et al, 1994], self-organization of gouge 
particles under large slip [Lockner and Byerlee, 1993; Scott et al, 1994; Scott, 1996], 
acoustic liquifaction [Melosh, 1996] and so on. In the second class of models, the fault 
has a low friction before the onset of the event. This may be due to the presence of 
low-strength clay minerals such as Montmorillonite (the weakest of the clay minerals) 
[Morrow et al, 1992], to an organized gouge structure similar to space filling bear- 
ings with compatible kinematic rotations [Herrmann et al, 1990], to the presence of 
phyllosilicates in well-oriented layers [Wintsch et al, 1995] or to the existence of a 
hierarchical gouge and fault structure leading to renormalized friction [Schmittbuhl et 
al, 1996]. Much attention has also been devoted to the role of overpressurized fluid 
[Lachenbruch, 1980; Byerlee, 1990; Rice, 1992; Blanpied et al, 1992; 1995; Sleep and 
Blanpied, 1992; Moore et al, 1996] close to the lithostatic pressure. The problem is 
that fluid pressure close to lithostatic value implies fluid trapping and absence of con- 
nectivity with larger reservoirs and the upper surface. Scholz [1992b] has also noticed 
that the fluid pressure scenario has one major problem that remains to be resolved, 
namely that a brittle material can never resist a pore pressure in excess of the least 
compressive stress a 3 without drainage occurring by hydrofracturing. Other sugges- 
tions to explain the absence of an anomalous heat flux near active faults include 
geometric complexity in the San Andreas fault at depth, hydrothermal circulation 
and missing energy sinks (see [Hickman, 1991] and references therein). 



3 Chemistry : mineral alteration and chemical 
transformation 

3.1 Alteration 

The first element of our approach consists in recognizing that the preparation for an 
earthquake starts at all scales, down to the molecular level. Under the action of a 
slowly increasing tectonic strain and in the presence of water which tends to concen- 
trate within fault zone discontinuities [O'Neil and Hanks, 1980; Thurber et al, 1997] 
and within defects in the minerals, rocks undergo a progressive hydration and lattice 
distorsion. Chemical reactions are slowly taking place within the rock fabrics. These 
processes involve the interaction of water molecules, both intact and decomposed 
into hydroxyl and hydrogen components by the exchange of electrons with the SiOl~ 
building block of silicates and also with the impurities within the rock minerals. For 
instance, phyllosilicates can precipitate in rock systems in the presence of water, 
provided magnesium is present [Wintsch et al, 1995]. Crystal plastic deformation 
(through dislocation movement) is enhanced by an atomic scale interaction of a com- 
ponent of the water with the Si — O bond structure. As a consequence, the minerals 
deform progressively, storing an increasing density of dislocations that are nucleated 
and stabilized by the presence of hydroxyl and other impurities in their cores. Some 
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deformation is displacive, i.e. corresponds to a distortion of the crystalline mesh as 
for instance in the local Si — O — Si unit flipping induced by the presence of a single 
water molecule [Jones et ai, 1992]. Other deformation is plastic and is due to the 
irreversible creation and motion of dislocations. These processes lead to a weakening 
of the rock materials. For small strains, this is nothing but hydrolytic weakening. 
We are interested in its extension in the large strain regime. For large strains, the 
effect of water is still largely unquantified and work on this problem is a priority 
in the search for a deeper understanding of earthquakes and more generally shock 
metamorphism [Nicolaysen and Ferguson, 1990]. 

This localized deformation and weakening of the rocks is preferentially concen- 
trated within discontinuities between rock fabrics since it is the domain where the 
largest amount of water is thought to be available for chemical and hydrolytic attack. 
Shear can then localize within a few ten-centimeter-thick ultracataclasite zone at the 
center of fault cores. These ultracataclasite zones are composed primarily of multiply 
reworked vein materials, testifying to repeating rupture and healing [Scholz, 1992]. 
In fact, preexisting faults are not necessary, only the preexisting heterogeneity in the 
mineral structures that can occur at many scales. The faults will appear as a con- 
sequence of the rupture on these discontinuities. Material discontinuities, textures, 
microcracks, faults are pervading the earth crust on all scales. It is important to 
recognize that physical fields (stress, strain, temperature, fluid pressure...) are, as 
a result, also highly heterogeneous. The formation of fault structures can occur by 
repeated earthquakes [Sornette et ai, 1994]. Earthquakes are localized on narrow 
zone in which the pressure, temperature and water content is particularly favorable 
for the preparatory stage. It is possible that faults are also partially selected and 
transformed by other processes than those occurring during earthquakes, for instance 
by ductile, creep and plastic localization [Ode, 1960; Orowan, I960]. 

Due to the weakening process and the plastic deformation and softening of the 
minerals, a large deformation may be progressively concentrated in a narrow domain 
while the intact parts of the rocks deform elastically at a smaller rate. This alteration 
does not need to occur in a constant steady state as it is controlled by the amount of 
available fluid which may be intermittent in time. Ultimately, alteration is coupled to 
the water content through a feedback loop involving its transport within the crust. If 
the localized deformation occurring at depth is coupled mechanically to the surface, 
it can be seen if an edifice or an object is sitting right across the discontinuity at the 
surface (see for instance [Bolt, 1993], pages 92-93) : the phases of observed continuous 
slips can be interpreted as the response of the upper gouge material to reactions in 
deeper rocks upon the introduction of water). One observes on some segments of 
the San Andreas fault important deformations concentrated in narrow meter-wide 
zones. On other segments or on other faults comprising about 90 % of all cases, such 
creep is not observed, probably due to the progressive decoupling (dampening) of the 
deformations by the overlayers. 

The total tectonic strain is not solely accomodated by a given localized zone. 
In addition to this localized deformation, large scale strain develops over the entire 
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loaded region due to tectonic motion. This large scale strain is the signature at 
large scale of the tectonic load on the elastic plates. The stress is transferred to 
the regions of local alteration and concentrated deformation through the agency of 
the large elastic domains which are thus necessarily deformed. It is thus important 
to understand that both styles of deformation coexist as a consequence of stress 
conservation. Only a fraction of the total deformation is accomodated in the core of 
the faults. 

3.2 Transformation of mechanical energy into chemical en- 
ergy during alteration via polymorphic phase transfor- 
mations 

Under the action of tectonic motion driven by thermal mantle convection, the crust 
is storing energy progressively. It also dissipates a part of the driving energy under 
the form of ductile and plastic deformations. If balance between input and output is 
positive as when the modes of ductile deformation are sub-dominent , more and more 
energy is being piled up in the crust. This cannot go on for ever and must be released 
suddenly in the form of an earthquake. In the broadest sense, an earthquake is the 
sudden release of a fraction of the stored energy, whatever the nature of this storage. 
In a nutshell, the earthquake cycle thus comprises the process of a slow energy storage 
ending in a brutal energy release. To understand what causes an earthquake, we must 
consider all possible significant forms of energy storage because the sudden event must 
release all the forms of energies that are not dissipated continuously. The form of the 
energy storage must in turn control the properties of the earthquake. 
In principle, the possible form of energy storage are 

1. elastic energy stored in the rocks as a result of their elastic deformation under 
the tectonic loading : this is the standard form considered in most models of 
earthquakes, as for instance in the elastic rebound model [Reid, 1910]. 

2. electric, dielectric or more generally electromagnetic energy : this form of energy 
storage is expected to be a very small fraction of the total stored energy due to 
the smallness of piezoelectric and electrokinetic conversion factors in the crust. 

3. "chemical" energy: we have seen that rocks are subjected to pressure, temper- 
ature, fluid and chemical conditions that put them in deformation regimes that 
are often much more complex than simple elasticity, especially close to rock 
discontinuities where fluid can have a very important role. In the literature, 
non-elastic modes of deformation are usually taken into account as factors that 
control the dynamics rather than the nature of end products. For instance, the 
nucleation and start-up of the earthquake rupture are often described in terms 
of non-elastic models of deformations. We propose here that non-elastic modes 
of deformations may play an important role also in determining the nature of an 
earthquake and not only its dynamics. We thus suggest that chemical energy 
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storage can be significant. Minerals in fault zones can be considered as analogs 
of piezo-chemical plants, which convert, in the presence of water and chemical 
environment, a part of the elastic energy into chemical energy that is stored in 
the form of new compounds. 

The proposed scenerio introduces the action of stress and strain in the stability 
of minerals. The proposed mechanism belongs to the class of processes known as 
mechanochemistry. Indeed, the motion of a dislocation by that of kink (leading 
to plastic deformation) is akin to a local chemical reaction in which an embedded 
"molecule" is dissociated, and then one of the product atoms joins with an atom 
from another dissociation to form a new "molecule" [Gilman, 1993]. Now, chemical 
reactions can be triggered by mechanical forces in solid phases, because unlike gases 
and fluids, solids support shear strain. Shear changes the symmetry of a molecule 
or a solid and is thus effective in stimulating reactions, much more so than isotropic 
compression [Gilman, 1996]. The reason for the strained minerals to be able to 
transform into metastable minerals lies in the kinetics. For application to the crust, 
we need to better understand the constraints in the parameter spaces (pressure, 
temperature, water affinity, impurities, strain, strain rate, rock composition) that 
may control the chemical transformations. 

Let us conclude this part by generalizing. There is growing recognition that 
mineral structures can form at much milder pressures and temperatures than their 
pure phase diagram would suggest, when in contact with water or in the presence of 
anisotropic strain and stress. For instance, diamonds can now be formed under rel- 
atively low pressure (100 MPa) and temperature (500 °C) under hydrothermal con- 
ditions [Zhao et al, 1997; DeVries, 1997], while the direct transformation route from 
carbon requires a pressure above 12000 MPa and a temperature of about 2000 °C. 
Another case in point is that the application of uniaxial stress along prefered direc- 
tions in quartz minerals results in the appearance of a new crystalline phase, where 
all silicon atoms are in fivefold coordination [Badro et al, 1996]. The stress threshold 
for the transition is lowered in this case by the application of the anisotropic stress. 
Novel behavior can also appear such as the instability of the melt-crystal interface 
[Grinfeld, 1986; Thiel et al, 1992] in the presence of nonhydrostatic strain. These 
recent discoveries are suggestive of the wealth of new phenomena that are possible 
when chemistry and/or phase transformations are coupled to anisotropic mechanical 
deformations. 

4 Dynamics : explosive release of chemical energy 
4.1 Energetics 

The metastable minerals are forming under tectonic stress and thus grow with miner- 
alogic orientations governed by those of the applied stress. They also start to deform 
until they become unstable and convert back to more stable minerals. We propose 
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that the back-reaction may become explosive, due the fast release of energy when the 
minerals become unstable due to the storage of dislocations and other defects. Only 
under these circumstances can an explosion occur. If the minerals are not severely de- 
formed, the reversion will be smooth as observed in some experiments [Green, 1972]. 
The experimental challenge, to test our ideas, is thus to produce sufficient deforma- 
tion after the nucleation of the metastable minerals to reach the metastable-unstable 
threshold at which a small disturbance explosively modifies its structure. A necessary 
condition is that the metastable phase should not recrystallize back spontaneously 
and progressively, as in coesite [Green, 1972]. 

Kuznetsov [1966] has determined the thermodynamic conditions for a transition 
from a metastable state of matter to a stable equilibrium state to occur as a detonation 
shock. He shows that many first-order (including polymorphic) phase transitions 
can lead in principle to a detonation. The general condition involves the position 
of the adiabat with respect to the isotherm in the pressure-volume diagram of the 
stable substance. Barton et al. [1971] and Hodder [1984] has taken up the idea and 
suggested its possible relevance for deep earthquakes. Randall [1966] has calculated 
the resulting seismic radiation from a sudden phase transition. 

To understand quantitatively and illustrate this process, we use the Landau theory 
of phase transitions in minerals [Salje, 1990; 1992; Heine et al, 1990]. Landau theory 
is simple and economic, as it uses symmetry constraints to derive the free energy 
as an expansion of the order parameters, thus minimizing the need for a detailed 
mineralogic description. Notwithstanding this simplicity, when the order parameters 
are correctly identified, it reproduces experimental observations of phase transitions 
in minerals with sufficient accuracy to be useful for many applications. It can also be 
extended to account for the dynamics of phase transitions [Salje, 1990; 1992]. This 
formalism should apply to any phase transition from one mineral form to another as 
occurs generally for instance in feldspars and other crystalline structures found in the 
crust. To keep the discussion as simple as possible, we use the simplest conceptual 
model of a coupling between a structural order parameter and the applied strain. In 
general, minerals and their structural phase transitions involve several coupled order 
parameters Qi, which describe the atomic displacements within the lattice structure. 
These parameters form a tensor that couples to the strain tensor e. 

Consider a structural transition described by the following free energy : 

G = -\aQ 2 + -bQ A + hQ + dQe + \gt 2 . (1) 

We consider the simplest case of a single scalar order parameter Q. More complex 
situations do not modify the mechanism and generalization of our discussion to more 
complicated situations is straightforward conceptually. The single order parameter 
Q represents for instance a pure dilational coupling or a pure shear coupling, ^ge 2 
is the elastic energy density of the material with elastic modulus g. The coefficient 
d quantifies the strength of the coupling between order parameter and strain. The 
parameters a and b are phenomenological coefficients for the phase transition. The 
"field" h controls the breaking of symmetry between the two phases. 
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In the absence of strain, the critical transition occurs at a = and separates a 
phase with Q = (for a < 0) from a phase with non-zero Q (for a > 0) possessing 
the Q — > — Q symmetry. For a > 0, h = and vanishing strain e = 0, the two phases 
Q± = ±-y/f have the same free energy. The symmetry between Q and — Q is broken 
by the field h. Keeping a > and varying h allows one to describe a first-order 
transition between the two phases Q + and Q_ with a jump in order parameter when 
h goes through zero. Note that there is no need of a field if we introduce a Q 3 term 
in the free energy expansion (|TJ) which breaks the Q — ► — Q symmetry. Other forms 
up to Q 6 have been considered to describe tricritical mineral phase transitions [Salje, 
1990; 1992; Heine et al, 1990]. Here, we use this simple expression as this is enough 
to demonstrate the effect. 

We now examine the influence of the coupling between the order parameter q 
and the strain field e. In fact, strain is recognized as an essential ingredient in 
structural phase transitions because elastic strain coupling is the dominant interaction 
between atoms and mineral cells in structural phase transitions [Marais et al., 1991; 
Salje, 1991; Bratkovsky et al, 1995]. For the sake of illustration, we make the stable 
undeformed mineral correspond to Q- and the metastable phase to Q + , which are 
obtained as the two minima of G. Note that the presence of strain e > has, in 
this model, the same effect as an increase of the symmetry-breaking field. Other 
forms of coupling do not have such a direct equivalence but exhibit however the same 
qualitative properties. In the presence of a finite strain e > 0, the two phases have 
their order parameter which are solution of ^ = (extremum of the free energy) 

with ^jf > (stability condition: local minimum). The first condition gives the 
cubic equation 

Q 3 + AQ + B = , (2) 
with A = f and B = *±^. Three cases occur [Beyer, 1991]. 

1. ^- + < : eq.(||) has three real roots, two of which are stable. This is the 
regime where the two phases Q- and Q + are locally stable and the phase Q + 
with the highest free energy is metastable. This condition holds for e < e* = 

_h j 2 a? 

d ~'~ y/Tf d ' 

2. ^- + Yf > i.e. e > e* : eq.(§) has one real solution and two conjugate 
complex solutions. Only the real one has a physical meaning and corresponds 
to the unique stable phase. 

3. ^ + = 0, i.e. e = e* : eq.(Q) has three real roots of which at least two are 
equal. For this critical value e* of the strain, the metastable state Q + becomes 
unstable as \q + vanishes and transforms into the stable phase Q_. 

Figure 1 schematically represents the whole process. Starting from a double well 
configuration with Q_ more stable than Q + , the deformation applied to the Q_ 
phase creates a higher energy state which eventually becomes comparable to Q + . 
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As a consequence, the system transforms into the metastable Q + . As the strain 
continues to increase, the free energy landscape deforms until a point where Q + 
becomes unstable and the mineral transform back into undeformed Q_. 

The energy released in the chemical transformation is estimated as follows. The 
free energy difference between the metastable phase and the stable mineral is taken of 
the typical order of 10 2 kJ /kg ~ 2 x 10 8 J/m 3 . Consider a volume of 40 km by 10 km 
by 0.1 m, corresponding to the core of a fault activated by an earthquake similar to 
Loma Prieta (1989 earthquake, M w = 6.9) [Kanamori and Satake, 1990] or Landers 
earthquake [Cohee and Beroza 1994], in which alteration has transformed about 1% 
of the minerals into the energetic metastable phase. This volume can then release 
10 14 J from the explosive phase transformation (compared to the energy released, 
say, by the Hiroshima bomb equal to 13 kT of TNT ps 5.2 x 10 13 J). For Loma 
Prieta, the seismic radiation energy is in the range [Kanamori et al, 1993; Houston, 
1990; Choy and Boatwnght, 1995] E s = 10 15 ~ 16 J. Another estimate is given by the 
static elastic calculation using a simple shear rupture model [Knopoff, 1958] giving 
a total released energy Et ~ ^Gd 2 L. This calculation assumes that the stress drop 
is equal to the average stress level prior to the event. For d ~ 1.6 m over a rupture 
length L m 40 km, we get E T p=j 1.3 x 10 15 . 

4.2 Explosive shock propagation 

The first chemical transformation is a slow process as it is fed by a slowly increasing 
deformation. In contrast, the transformation from the deformed metastable phase 
into underformed minerals is not slow and has a dynamics which is expected to be 
linear in time (i.e. a front velocity can be defined). The reason is the following. 
If a structural phase transition occurs at non-constant chemical composition, the 
dynamics is in general diffusive and thus slow. However, if the composition is constant 
(as is expected here for structural phase transitions in minerals), the transformation 
may occur as a front propagation (so-called "massive" transformation) [Christian, 
1965]. Among massive transformations, the displacive transformations can even occur 
in volume at an extremely fast rate since it involves only local bond rotations (as in 
quartz a — > quartz (5 for instance). This last situation does not apply in general to 
structural phase transformations that require rupture of bonds and not only atomic 
rotations. A "massive" transformation is obtained by the propagation of curved, 
flexible boundaries which move with variable speeds and have the ability to cross 
grain boundaries. The velocity is orders of magnitude higher than that of a reaction 
involving long-range diffusion. During transformation, migration over only a few 
interatomic distances is required. 

There is a more fundamental reason for the fast propagation of the explosive front. 
The point is that the products left in the wake of the front are highly fragmented 
minerals and not the well-structured arrangements that would be the result of a slow 
diffusive-limited phase transformation. The usual slow velocity of phase transfor- 
mations is due to this diffusion-limited nature which is absent here in an explosive 
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process. 



If the transformation can be slow or fast depending on the pressure and tem- 
perature conditions, it has been shown that the chemical reaction at the atomic 
level is not essentially different even when explosive [Lonsdale, 1969]. In the mineral 
transformations proposed here, the dynamics should be even faster than for stan- 
dard "massive" transformations in alloys. The reason is that no thermal activation is 
needed as the free energy barrier is made to vanish by the increase of strain. At the 
atomic level, large strains allow a derealization of electronic charges helped by impu- 
rities which lead to a new mineral structure [Gil-man, 1992; 1993; 1995b; 1996]. This 
is a situation which is very similar to what happens with solid explosives. Indeed, 
explosive substances are nothing but species storing chemical energy in metastable 
chemical configurations that is released suddenly. Here, an energetic solid substance 
(the metastable phases) releases energy quickly by transforming to a more stable low 
energy substance. It has been proposed [Gilman, 1995a] that intense deformation 
by bending of atomic bonds occur in a very narrow zone of atomic scale that can 
propagate at velocities comparable to or even higher than the velocity of sound in 
the initial material. 

Let us show that the propagation of the explosion front occurs at supersonic 
speed. Let us consider an metastable phase that has grown aligned with a coherent 
crystalline structure under the influence of the global stress and strain fields. Imagine 
the following simple model. A one-dimensional chain is made of atoms of mass m 
linked to each other by energetic links of spring constant k which, when stressed 
beyond a limit, rupture by releasing a burst of energy Ag converted into kinetic 
energy transmitted to the atoms. Initially the chain of atoms is immobile. Suppose 
that the first atom on the left is suddenly brought to a position that entails the 
rupture of the first bond. This rupture releases the energy Ag that is converted 
into kinetic energies of the atom fragment that is expelled to the left and of the next 
atom to the right which becomes the new left-boundary of the chain. All atoms along 
the chain start to move progressively due to the transmission of the motion by the 
springs. Now, due to the impulsive boost y / ^f that the boundary atom received, 
it will eventually stress the bond linking it to the next atom towards its rupture 
threshold. When this occurs, it is expelled by the energy that is released and the 
next atom forming the new boundary is itself boosted suddenly by the amount J— . 




It is then clear that this leads to a shock propagating at a velocity larger than the 
sound velocity since the atoms are receiving boosts that accelerate their motion 
faster than what would be the propagation by the springs with the usual acoustic 
wave velocity. Taking the continuous limit, the resulting supersonic shock velocity U 
is given by an adaptation of the formula proposed by Gilman [1995] : 



where c is the longitudinal (P-wave) velocity (around 5000 m/s), and AG ~ 
10 2 kJ/kg is a typical value for the free energy release by the transformation. This 
yields U « 1.003 c. 




U 2 = c 2 + AG 



(3) 
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We present in the appendix 1 a calculation of the properties of such a front prop- 
agation during an explosion in a one-dimensional tube configuration. The calculation 
highlights the dependence of the front velocity on the energy released by the explo- 
sion. The appendix 1 treats a solid-solid explosion with a front separating a material 
with high energy ahead of it from a region of denser and strengthened material with 
lower energy behind it. The front propagation is found supersonic with respect to the 
wave velocity (here in ID) in the weaker material ahead of it but remains subsonic 
with respect to the wave velocity behind it. 

Associated to the mechanical transformation and the associated elastic waves, we 
also expect electrical signals to be generated. Indeed, the metastable phase destabi- 
lization and transformation to other minerals occurs, as already stressed, in coherent 
aligned crystals. As a consequence, a net non-vanishing piezoelectric effect should 
appear and the transformation to other minerals must induce a significant electro- 
magnetic pulse. The elasto-electric coupled modes in a piezoelectric material are 
calculated in the appendix 2. The appendix 3 revisits the calculation of the appendix 
1 in the presence of the electric coupling process. 

5 Dynamics : the genuine rupture 

5.1 Explosive fluidization and unlocking of the fault 

The explosive transformation of the deformed unstable energetic mineral phase is 
violent and leads to fragmentation with the generation of intense shaking due to 
high frequency sound waves that remain trapped in the loosened low acoustic impe- 
dence core of the fault. We argue that this leads to acoustic fluidization [Russo et 
ai, 1995]. Melosh [1996] has proposed acoustic fluidization as a mechanism for the 
low dynamical friction of faults (the initial unlocking of the fault is not described in 
his scenario). We differ from him in that the source of the high frequency acoustic 
modes is not the rupture propagation but rather the explosive chemical phase trans- 
formation. Furthermore, the acoustic pressure does not need to reach the overburden 
pressure to produce a fluidization of the fault. Sornette [1998b] has shown that there 
is a problem with Melosh's mechanism because it predicts a slip velocity during an 
earthquake more than two orders of magnitudes smaller than the typical meters per 
second for observed earthquakes. 

It is possible to save this mechanism by invoking that the acoustic pressure does 
not need to reach the overburden pressure but only a small fraction 77 of it, in order to 
liquidity the fault. Indeed, it is well-established experimentally [Biarez and Hicher, 
1994] that the elastic modulii of granular media under large cyclic deformations are 
much lower than their static values. This effect occurs only for sufficiently large 
amplitudes of the cyclic deformation, typically for strains e a above 1CT 4 . At e a = 1CT 3 , 
the elastic modulii are halved and at e a = 10 -2 , the elastic modulii are more than five 
times smaller than their static values. As a consequence, the strength of the granular 
medium is decreased in proportion. Extrapolating these properties to the crust, we 
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need to estimate the strain created by the acoustic field. The acoustic pressure is 
related to the acoustic particle velocity v by p = pcv. Assuming p = rjpgh, this 
yields v = ngh/c pa 12 m/s for p pa 200 MPa, a density p = 3 10 3 kg/m 3 , a velocity 
c = 4000 m/s and rj = 0.1. At a frequency /, this corresponds to an acoustic wave 
displacement u a = v/2irf Pi 2 10 -3 m at / ~ 10 3 f/z. The corresponding strain 
u a /w is ~ 2 10~ 3 for a gouge width u> of the order of one meter [Melosh, 1996] over 
which the intense shaking occurs. These estimations suggest that Melosh's criterion 
that the acoustic stress fluctuations must approach the overburden stress on the 
fault for acoustic fluidization to occur is too drastic and smaller shaking can reduce 
significantly the fault friction. 

The peak deformation strain amplitude e a of the acoustic waves generated by the 
explosion can be estimated from the density change emsMstaMe^pstabie ^hat we ^ a ke Q f 

Pstabie 

the order of 10% during the the explosive transformation. This leads to e a pa 3 — 4 %. 
The acoustic pressure equals to e a pgh, i.e. a few percent of the lithostatic pressure. 
This is almost two order of magnitude less than the value needed in the acoustic 
fluidization mechanism proposed by Melosh [1996]. In the presence of this shaking, 
the strength of the granular gouge is drastically decreased and the gouge can slip 
under the applied tectonic stress. In addition, if the stress-strain curve of the gouge 
has a maximum at the pressure and temperature conditions of the seismogenic depth, 
usually occurring at a deformation of a few percent, the acoustic waves created by the 
explosive transformation may even lead to an intrinsic shear instability, analogous to 
the localization instability in sand. A third fluidization mechanism may also add to 
the instability. When the explosive transformation occurs, the ensuing shaking starts 
to deform the gouge material. The first response to deformation of granular material 
is to compress, even if it becomes dilatant at larger deformations. As a consequence, 
any interstitial fluid is first compressed which decreases the friction force between the 
grains and thus may lead to a release of the fault, which in turn may start to rupture. 
This effect is known as the (standard) fluidization of granular media. 

5.2 Rupture propagation and seismic radiation 

There are two possible contributions to seismic radiation. First, the phase transfor- 
mation by itself radiates seismic waves: Randall [1966] and Knopoff and Randall 
[1970] have shown that a sudden change of shear modulus of the material in the pres- 
ence of a pre-existing shear strain leads to a "double-couple source" with the correct 
and usual characteristic radiation pattern found in earthquakes at low frequencies. In 
fact, it is impossible to discriminate between this source and that of the displacement 
dislocation from their first motions. The size of the corresponding double-couple 
moment is [Knopoff and Randall, 1970] 

M — 28 p: 5e V , (4) 

where 5p is the change in shear modulus, Se is the change in strain inside the volume 
V in which the phase transformation has occurred. As an order of magnitude, we 
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take a change in modulus of the order of 10 10 Pa and an upperbound for the change 
in strain 5e one-third the relative change in density (10 %). The volume V is the 
surface S of the fault times the width w ~ 0.1 meter of its core over which the phase 
transformation occurs. We thus see that this corresponds to the seismic moment of 
an earthquake with average slip w 5e equal to a fraction of centimeter, irrespective of 
the size of the fault. This contribution is thus negligible for large earthquakes having 
slips of meters or more, but may become significant for small earthquakes which are 
also well-recorded. In addition to this double-couple component, the change in bulk 
modulus and in density radiates isotropically with radial directions for the first P 
motions and no first S motions. However, this is again a small effect. 

The second contribution to seismic radiation comes from the mechanical fault 
slip. The shaking of the fault core induced by the explosive transformation unlocks it; 
indeed, the fluidization implies that the fault can no longer support the initial loading 
stress, since its strength tends to vanish. This explosive transformation triggers the 
fault slip and the fault starts to slip under the action of the pre-existing shear stress 
and radiates seismic waves. A crucial point in this model is that the fault slip is not 
triggered by reaching a stress threshold (corresponding either to friction unlocking or 
rupture nucleation), but rather by a chemical instability : as a consequence, any level 
of stress will activate the fault slip when the explosive phase transformation occurs. 
The explosion could also be triggered from another earthquake. This possibility needs 
further investigation. 

The rupture dynamics is controlled by the usual elasto-dynamic equations. The 
rupture propagation lasts as long as the high-frequency waves that are trapped within 
the shaking gouge with low-acoustic impedence [Harris and Day, 1997; Li et al, 
1994] remain of sufficient amplitude to unlock the fault. The detrapping of these 
waves control the healing of the fault, and therefore the static stress drop. This 
mechanism is expected to produce large variations of static stress drops, depending 
on the degree of chemical alteration and storage of chemical energy necessary to 
obtain the unlocking of the fault. This is controlled by the detailed mineralogy in the 
fault core and the availability of fluids. The total size of the rupture is determined by 
the extension of the domain over which the supersonic shock has propagated, which 
is itself controlled by the alteration processes that have matured the material and 
stored a suitable amount of chemical energy. A fault region that is weakly or not 
altered plays the role of an energy sink for the explosive shock propagation and will 
tend to stop it. 

6 Consequences and predictions 

The proposed theory suggests to explain a number of observations into a coherent 
framework. 

• Strain, stress and heat flow paradoxes : There is no need for elastic strain con- 
centration over a scale of about 10 km (which, as we have reviewed, is usually not 
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observed) and very localized plastic-ductile strains are expected. There is no need for 
large stress to unlock the fault and the low friction is generated dynamically, prevent- 
ing heat generation and providing a solution to the heat flow paradox [Lachenbruch 
and Sass, 1980]. 

• The longer the recurrence time, the larger the stress drop : Several studies have 
shown recently a remarkable relationship between the average slip rate on faults and 
the stress drop associated to earthquakes occurring on these faults (see [Kanamori, 
1994] for a review). Earthquakes on faults with long repeat times (thousands of years) 
radiate more energy per unit fault length and have a significant larger dynamical stress 
drop than those with short repeat times (a few decades to a few centuries). In our 
framework, a longer period gives more time to saturate the chemical energy storage 
and leads to more "energetic" earthquakes. 

• Seismic P-wave precursors: Seismic P-wave "nucleation phases" have been re- 
ported [Beroza and Ellworth, 1996] that seem to precede the arrival of the first P-wave 
(longitudinal compressive acoustic wave). These observations are still controversial 
[Mori and Kanamori, 1996], not only because the reported signals are weak and 
the effect is hard to establish, but also because their presence is essentially ruled out 
within the standard pictures. The proposed explosive mechanism and resulting shock 
wave propagating at slightly supersonic velocity (1.003 c according to our estimate) is 
a natural candidate to rationalize the observation of these seismic P-wave precursors, 
if they exist. We predict an advance of about 6 milliseconds of the precursor to the 
first P-wave motion for a propagation over 10 km between source and a seismic sta- 
tion. This seems of the correct order compared to observations [Beroza and Ellworth, 
1996]. However, we predict a delay proportional to the distance to the station while 
Beroza and Ellworth [1996] find a delay proportional to the fault rupture length. The 
finding that the duration of a precursor scales with the size of the earthquake is in 
agreement with our model in which the mechanical rupture occurs on the length of 
the fault over which the explosive shock has occurred and thus we should expect the 
precursor to be proportional to the total size of the rupture. 

• Tilt anomalies have sometimes been reported before earthquakes (for instance 
before the Haicheng earthquake [Scholz, 1977; Kanamori, 1996]). The magnitude of 
the tilt is usually very small (1 — 20 fir ad, say) and not always present. The slow 
phase transformation during the alteration process and chemical energy storage leads 
to a density change that produces a weak surface deformation whose direction and 
amplitude depend on the mechanical heterogeneity of the earth. This may suggest 
a source for the tilt that is sometimes observed. For a phase transformation where 
10% of the minerals undergo a relative density change of 10 % occurring within the 
fault core over a depth of 10 km and a width of 1 meter, the expected tilt anomaly 
is ~ 3 firad over a distance of 10 km. 

• Earth tides: The long-term tectonic loading stress rate of the order of 
10~ 3 bar /hr is much less than the stress rate up to 0.15 bar /hr due to earth tides from 
gravitational interactions with the Moon and Sun. Tidal triggering of earthquakes 
would thus be expected if rupture began soon after the achievement of a critical 
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stress level. The most careful statistical studies have found no evidence of trigger- 
ing [Heaton, 1982; Rydelek et al, 1992; Tsuruoka et al, 1995; Vidale et al, 1998]. 
Two existing theories can explain these observations by invoking high stress rates 
just before failure. Dieterich's model of state- and rate-dependent friction predicts 
high stressing rates accross earthquake nucleation zones [Dieterich, 1992]. Alterna- 
tively, changes in fluid plumbing of the fault system could conceivably be more rapid 
than tidal strains and may trigger failure [Sibson, 1973]. Our theory, which does 
not attribute the triggering of an earthquake to a critical stress threshold, is fully 
compatible with these observations. 

• Seismicity remotely triggered at long distances: There is strong evidence that 
the Landers earthquake, June 28, 1992, southern California, triggered seismicity at 
distances equal to many times the source size [Hill et al, 1993], with a rate which 
was maximum immediately after passage of the exciting seismic waves. The problem 
is that the dynamical stress created by the seismic waves is very small at these long 
distances, of the order of the effect of lunar tides (0.01 MPa and less) which have not 
been found to be correlated with earthquakes. Sturtevant et al. [1996] have recently 
proposed a model in which the earthquakes are triggered by a rapid increase of pore 
pressure due to rectified diffusion of small pre-existing gas bubbles in faults embedded 
in hydrothermal systems. This model depends on the confluence of several favorable 
conditions, in particular supersaturated gas in water, large mode conversion occurring 
in the geothermal field and a very low (50 m/s) shear velocity in the porous media 
filled with the bubbly liquid. Alternatively, the chemical instability we propose should 
be much more sensitive to high frequency waves than to low frequency modulations 
(think of the jerky motions that one strives to avoid when manipulating nitroglycerin 
explosives!). 

• Pre-seismic chemical anomalies : Various anomalous precursory chemical emis- 
sions have been reported recurrently, but not systematically as for all other proposed 
precursory phenomena. Nevertheless, it seems to be a phenomenon worthy of study. 
For instance, Sato et al. [1986] reported anomalous hydrogen concentration change 
in some active faults, and in particular in association with the Coalinga earthquakes. 
Tsunogai and Wakita [1995] reported anomalous ion concentrations of groundwater 
issuing from deep wells located near the epicenter of the recent earthquake of mag- 
nitude 6.9 near Kobe, Japan, on January 17, 1995. Similar anomalies have also been 
measured for Radon emission [Igarashi et al, 1995] for the same earthquake. Jo- 
hansen et al. [1996] presented a thorough statistical analysis of these precursors for 
the Kobe earthquakes and concluded that these time- dependent anomalies are well 
fitted by log-periodic modulations around a leading power law [Sornette, 1998a], qual- 
ifying a critical cooperative behavior. The source of these chemical anomalies may be 
different but seem to point to a chemical source, that we tentatively associate with 
the chemical transformations that may accelerate close to the instability. 

• Electric effects: The metastable crystals are expected to form with a prefered 
coherent orientation. For those minerals that do not present centers of symmetry, 
they should thus exhibit a net piezoelectric effect. Their explosive transformation 
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may then lead to significant co-seismic electric signals. Precursory electric signals 
might be associated with partial precursory conversions of the metastable minerals. 
This might provide a scenario for rationalizing precursory observations [Debate on 
VAN, 1996] and help in improving their investigation. 

• Deep earthquakes: Finally, our theory suggests that deep earthquake, which 
have been proposed to be due to unstable olivine-spinel transformations [Green and 
Houston, 1995], are not so different from superficial earthquakes, not only in the 
nature of their seismic radiation but also in their source mechanism. This may call 
for a re-examination of the phase transition mechanism for deep earthquakes in the 
light of the action of water and other impurities in presence of finite strain. 

• Inversion of metamorphic data : The usual inference of past tectonic conditions 
from the examination of minerals and their corresponding equilibrium thermody- 
namic phase diagrams may be questioned in view of the evidence summarized here 
on out-of-equilibrium processes in the presence of water and finite strain. This may 
call for a reexamination of the models of crustal motions based on inverting meta- 
morphic patterns in fault zones. 
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7 APPENDIX 1 : Explosive shock neglecting elec- 
tric effects 

Our treatment follows the analysis of Courant and Friedrichs [1985] of shock waves 
and we adapt it to the case of solid phase transformations. We restrict the treatment 
to a one-dimensional system. In this simple representation, it does not make difference 
in the mathematical description whether we consider a compressive, extensive or shear 
(antiplane deformation). The formalism below ressembles the treatment of Courant 
and Friedrichs [1985] for wave propagation of finite amplitude waves in elastic-plastic 
materials, but differ in one essential point, namely in [Courant and Friedrichs, 1985] 
the stress-strain characteristics of plastic material is of the weakening type and thus 
does not allow for shocks. We study the opposite case where a shock occurs. 

We consider a long bar parallel to the Ox axis. The bar is initially deformed 
uniformly with a strain s T which is the sum of an elastic part si and a ductile- 
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plastic part Si : s T = si + sf. In the first initial solid phase, the elastic modulus 
is G\. A perturbation is brought at one extremity of the bar and a new crystalline 
phase is nucleated and the phase transformation propagates along the bar. The new 
solid phase has a different larger elastic strain s 2 and thus smaller plastic strain 
sf = s T — S2, since the total strain remains fixed. The structural transformation can 
thus be viewed as a transmutation of plastic into elastic strain by the local atom 
rearrangements within the crystalline mesh. The new phase has a different stronger 
elastic modulus G2- This transformation is thermodynamically favorable because the 
total energy, sum of the internal energy and of the elastic energy, decreases even if 
the elastic energy increases. The plastic part is taken into account by the internal 
thermodynamic energy which determines what is the stable solid phase. 

Our goal is to calculate the characteristics of the propagation of the phase trans- 
formation which will turn out to be a shock. For the sake of being specific, we consider 
an extensive strain (the conclusions are identical for compressive or shear strains). 
We will be concerned only with the elastic part of the strain as it is the only one 
which contributes to the elastic energy. This is done to simplify the treatment which 
is presented as a plausibility demonstration. This assumption amounts to neglect 
the variation of the fraction of the density accomodated by the plastic part of the 
deformation. All our presentation below thus substracts the plastic deformation. 

If a is the initial position of an atom in the bar, it becomes x(a, t) under some 
elastic deformation. The elastic component of the strain is given by 

If po is the initial density, conservation of mass reads p^da = pdx and thus 

— = — = 1 + s . 6 

p da 

The velocity of an atom along the bar is 



The equation of motion is 



du da 



^ dt dx 



The r.h.s. of (H) can be written lf = ifir = ¥¥¥- It makes sense to consider the 

v Ox aa ox as aa ox 

derivative of the stress with respect to strain, as we consider only the elastic strain 
and impose fixed plastic deformations in the two solid phases. Using the relation 
ffi = 1 + s and (|^), we get the wave equation 

d 2 x 2 d 2 x , . 

W = 9 da^ ' [ } 
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with 

/ 1 da . , 

» s v*& ■ (10) 

Note that g, c alled the rate of change or shift rate, is usually different from the sound 
velocity v = ^ — It is straightforward to check that 

Po9 = pv , (11) 

which implies that g increases when the material is denser. In the elastic range 
where @f = G is constant, we get g = ,/^- which is constant and v = ^-g which 



ds ' V Po P ' 

varies with the density p and thus with the deformation. Note that in the limit of 
small deformations, these differences can be neglected and g and v become indentical. 
Here we keep the distinction as it is important for finite deformations as occurs in 
"massive" structural transition in which the atoms move over a distance equal to a 
finite fraction of the lattice mesh. 

The structural solid-solid transition is modelled in this framework through the 
form of the function cr(s). Before the transformation, we assume a relatively weak 
elastic solid a = G\s. In the neighborhood of a deformation threshold s*, the modulus 
crosses over to a larger one G 2 and the characteristics becomes a = G 2 s for s > 
s*. This schematic dependence summarizes the nature of the solid-solid transition 
between a weak solid with elastic strain s\ to a stronger solid with larger strain s 2 . 

Consider a bar of material deformed uniformly with strain sq everywhere along Ox. 
Suppose that a localized perturbation or inhomogeneity produces a local deformation 
larger than s* at the left boundary of the bar. This perturbation is taken to represent 
the local nucleation of the stronger solid phase. The question we address is that of 
the growth of this new phase. Qualitatively, the density perturbation will start to 
advance to the right in the lighter phase. Since the phase rate g(s) increases with 
s, the largest deformations propagate the fastest. An initial smooth disturbance will 
progressively steepen and a shock will eventually form. The essential condition for the 
formation of a shock is thus the increase of g(s) with s. We imagine the extremity 
of the bar to be suddenly extended and to move with constant velocity Up in the 
immobile crystal ahead. No matter how small up is, the resulting motion cannot be 
continuous because a continuous motion would imply a forward-facing simple wave, 
i.e. a centered simple wave, in order to achieve a discontinuous change of velocity at 
the origin. However, the material velocity through a centered simple wave becomes 
negative if it vanishes ahead of the simple wave. Therefore, no adjustment to the 
positive piston velocity is possible by continous motion. The answer to this problem 
is that a shock front appears, moving away from the extremity of the bar with a 
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constant supersonic speed (with respect to vi = y uniquely determined by the 
density and the velocity of the quiet crystal and the piston speed. We note that 
the shock condition is not obeyed in the analysis of Courant and Friedrichs [1985] 
of wave propagation of finite amplitude waves in elastic-plastic materials. The shock 
studied here is indeed due to the rigidifying condition which is absent in elastic-plastic 
materials. 
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In the presence of a shock discontinuity, we have to write the conservation equa- 
tions in integral form as the usual differential formulation is not adapted to treat 
the discontinuity at the shock. We note ai(t) (resp. a 2 {t)) a point on the left (resp. 
right) side of the shock. 

The conservation of mass reads 

d r a 2{t) 

— pdx = . (12) 

dt Jai(t) 

The conservation of momentum reads 
d r a 2(t) 

— / pudx = p(ai,t) — p(a 2 ,t) = a(a 2 ,t) — a(ai,t) . (13) 

dt Jai(t) 

p is the external imposed pressure, which is opposed in sign by the internal stress. 
The conservation of energy reads 

d rai{t) ]_ \ ug 

— P(e + -w+- — )dx = p(a l ,t)u(a 1 ,t) -p(a 2 ,t)u(a 2 ,t) , (14) 
dt Jai(t) 2 2 p 

where e is the internal energy of the crystal. The r.h.s. corresponds to the work of 
the external pressure at the extremities of the bar. 

If all fields are continuous, we retrieve the usual equations of motion. Here, we 
assume, from the preceeding considerations, that a point of discontinuity exists within 
the bar at position x — £(t) and we note the velocity of the front 

f t =U(t) . (15) 

Consider an integral of the type J = j^l ty(x, t)dx, where $ is discontinuous at 
x — £ . Then, as shown in [Courant and Friedrichs, 1985], ^ must be calculated as 

dJ d /■£(*) d /•°2(t) 

M=itL nx ' t)dx+ dtL * (x ' t)ix > (16) 

with keeping in mind that ^> is discontinuous at x — £. In the limit where a\ — > a 2 
while still keeping the condition a\ < £ < a 2 , we get the fundamental shock equation 

^ = m 2 V 2 - VxVz , (17) 
dt 

where \J/i( 2 ) = ^(ait 2 ), t) and 

V m = u m - U . (18) 

Applying ( |TTD to the above conservation equations leads to the usual fundamental 
shock equations : 

PiVi = P2V2 = m , (19) 
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where m is the mass flux across the shock. 



mu\ +px — mu 2 + p 2 = c , (20) 

where c is the impulse. The equation of energy conservation yields, after some ma- 
nipulations, 

ei + 7 ;V 1 2 + ^ + -— = e 2 + -V 2 2 + ^ + -^ . 21 
2 pi 2 pi 2 p 2 2 p 2 



Using pO) and (pi), we find that 



™ 2 = -^? , (22) 



Pi P2 



which reduces to the wave velocity in the limit of very weak shocks. We thus recover 
the fact that a sound wave can be interpreted as an infinitely weak shock. This must 
in fact be obvious since the conservation equations then recover the wave equation 
directly. 

Using p = —a, the relation (|6]) between the strain s and density p and the two 
linear elastic Hooke's law (ai = Gis and a 2 = G 2 s) in the two phases on each side of 
the shock, we rewrite the energy equation (|2T| ) as 



e (p 2 ) + ^(4-3^ + 2) = e(p l) + ^(4-3^ + 2) . (23) 
2p 2 p 2 p2 2pi pf pi 



The two different solid phase structures of the crystal give us e(pi) and e(p 2 ) (or their 
difference). We can usually also determine the elastic coefficients G\ and G 2 of the 
two mineral phases. The equation ( p3|) thus determines p 2 as a function of p\. From 
this, we deduce the flux m from (|22|). And with (0), we get V\ and V 2 and then 
deduce u 2 (assuming u± = corresponding to a material at rest ahead of the shock) 
from (p0|). We then deduce the shock velocity U = — V\ = u 2 — V 2 . 

Under extension and with our conditions for the shock that < si < s 2 and 
o~i < a 2 , we verify that m 2 given by ( p2|) is positive. Indeed, m 2 can be written m? = 



Po °" 10 " 2 > for compression or extension as the sign results from the strengthening 



character of the transition, i.e. ^ > 0. Note that a usual ductile rheology has the 
opposite sign and therefore cannot develop a shock [Courant and Friedrichs, 1985]. 

Using ( pT8| ) together with (|I9"D, and with the initial repose condition that u± = 0, 
we get the formula for the shock velocity 



u= H = na_ 1 -a 1 
Po V Po s 2 - Si 

where we take the absolute value as m is negative, since the flux of matter goes 
from the right to left. We can compare the shock velocity U with the sound wave 
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velocities V\ and v 2 given by vi^) = \J~ %1 < '^ ) = "pcf"' ^ * s eas i er to compare 

p o jj2 _ g 2 s2-Gisi w ^ p oV z = Gi(2)(si( 2 ) + l) 2 . We construct the difference 

p U 2 - p v 2 1{2) = 

(_1)2(1) , 2 3 2 

- _ - ( s 2(i) [G 2 (i) - G 1{2 ) 0i( 2 ) + l) 2 ] + Gi( 2 ) (s? (2 ) + 2s 2 (2) 

We verify that p U 2 —p v 2 > for a very large set of parameters while PqU 2 — p§v\ < 
is realized if and/or ^ are sufficiently large. This condition is realized for instance 
with G 2 = 2Gi and s 2 = 3si = 0.3 while if S2 = 2s± = 0.2, it is not. The shock 
regime corresponds to the situation where v\ < U < v 2 . 



8 APPENDIX 2 : Elastic-electric coupled wave 

For the simplicity of the exposition, we restrict to a ID model and express the elastic 
and piezoelectric equation in scalar form. This does not preclude from the fact that 
piezoelectricity occurs in materials lacking a center of symmetry. Generalization 
to the full tensorial expressions involves straightforward manipulations. The elastic 
material obeys Hooke's law relating the stress a to the strain S : 

a = GS , (25) 

where G is the elastic modulus. The electric analog relates the electric induction D 
to the electric field E : 

D = eE , (26) 

where e is the dielectric coefficient of the material. In the presence of a non-vanishing 
piezoelectric effect, (p5|) and (|2~6"D become coupled through the piezoelectric equations : 



S = G- l a + dE , (27) 

D = da + eE , (28) 

where d is the piezoelectric coupling coefficient (= 2 10~ 12 CN^ 1 for quartz). 

In order to derive the propagative modes, we supplement these two constitutive 
equations ( P7| , P8"D with the fundamental elastic and electric equations. 



p— = — (29) 



piezoelectric coupling and using S = lr, we get the standard wave equation 



where u is the displacement at position x. Together with (^5|) valid in the absence of 

du 
dxi 

ld 2 u d 2 u m 



v 2 dt 2 dx 2 
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where the acoustic wave velocity is 



" S V7 • (31) 

The two Maxwell equations are 

dB 

rotE = -— , (32) 
1 - - dD 

—votB = — , 33 
fi at 



where \i is the magnetic permittivity of the material. Together with fl2l]| ) valid in the 
absence of piezoelectric coupling, the elimination of B yields the wave equation 



1 d 2 E 



c 



2 



dt 2 



AE , (34) 



where the electromagnetic wave velocity is 

1 (35) 

In the presence of the piezoelectric coupling given by (|27| , |28|) , we get the two coupled 
equations (keeping the single non-vanishing electric field component) 

d 2 E 1 M G j2 ,d 2 E d 3 u 

-{l--d 2 )— + dGii—— , (36) 



dx 2 c 2 e dt 2 dxdt 2 

d 2 u _ }_<?E_ £E_ 

dx 2 v 2 dt 2 dx 

We recover (|30| ) and (|34] ) in absence of coupling d = 0. The strength of the elastic- 
electric coupling (EEC) wave is measured by the dimensionless parameter — d 2 . For 
quartz, e = 4.5e and we take the value G w 30GPa. With d = 2 10" 12 CA^ _1 , we 
get — d 2 w 3 10~ 3 , i.e. a small coupling. 

We look for propagative modes E = i? e lfcx_w *, u = M e lfcrr_w *. Inserting in (|37| , |36|) , 
we obtain the dispersion relation, by the condition of vanishing the determinant, 
under the form of a quadratic equation in uo 2 as a function of k 2 : 

(l-^V-( C 2 + ^ )A; 2 a; 2 + c 2 u 2 A; 4 = (3g) 

e 

whose general solution is 



(c 2 + v 2 )k 2 ± V A 
2(1- ^d?) 



4 = V n/1 7 G * , (39) 
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where A = (c 2 + v 2 )k 4 — Ac 2 v 2 k 4 (l — ^d 2 ). Since the coupling ^-d 2 « 1, we can write 
with a good approximation A = k 4 (c 2 — v 2 ) 2 . Within this approximation, there is a 
fast mode propagating at close to the velocity of light and a slow mode propagation 

2 

slightly slower than sound. To first order in ^ , we get 

= c 2 k 2 \ + ^J « c 2 k 2 (l + -d 2 ) , (40) 
1 — —d l e 

e 

and 

^ = v 2 k 2 l + ^ a ^ ^ v 2 k 2 {l + (41) 

1 — —a 2 e 

Note that is real since we have neglected dissipation. From (|36| , |37|) , we get the 
relationship between the electric field amplitude and elastic displacement of the two 
modes : 

vi = -^\< • («> 

^ Gdk . 
£- = . (43) 

For a similar elastic amplitude Uq —Uq, we see that ^ = ^-^j << 1. This shows 

that the fast mode is essentially electromagnetic with negligible elastic deformation 
and the slow mode is mainly elastic with very weak electromagnetic fields. 

An arbitrary perturbation of the displacive type will decompose onto these two 
modes which will thus be excited and propagate in the narrow zone where minerals 
have been aligned to create a net piezoelectric effect. In particular, the amplitude 
of the electromagnetic mode is given by the equation (|28|), namely eE = da (in 
absence of preexisting electric field). Using e = 4.5eo with eo = 8.85 10~ 12 Fm^ 1 
and d = 2 10- 12 CJV -1 , we obtain E ~ 0.05a (Pa). For a stress of lOMPa, this 
yields an electric field of 5 10 5 V/m. This numerical value is of course only an 
order of magnitude as several uncertainties control its determination, such as the 
piezoelectric coupling d which is probably significantly over-evaluated and the stress 
amplitude. Nonetheless, this calculation suggests that significant electric fields can 
be created locally. 



9 APPENDIX 3 : Structural shock including elec- 
tric effects 

It is the nature of the structural transition of polarized atoms within the crystalline 
structure that electrons are released and exchanged during the propagation of the 
shock. We can even infer that, since the hydrolytic weakening is fundamentally a 
redox process, the process of the exchange of electrons is a key participant in the 
solid-solid phase transition. The release of electrons acts as a self-fuelling process by 
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the polarization source and local crystal structure distorsion it induces. Associated 
to the elastic shock wave, the modification of the crystalline structure induces a 
polarization wave, since electrons are released and exchanged. 

The mathematical treatment of the solid-solid phase transformation in terms of a 
shock presented in the appendix 1 can be extended to take into account the electro- 
magnetic coupling. The equations (|19|) of mass conservation and ( p0|) of momentum 
conservation remain the same. The equation ( pTj) of energy conservation is modified 
into 

1 T ^9 Pi 1 ClSl 1 -D? l Tr n Vo 1 CToSo 1 -Do , , ,, 

ei + ^ + - + ^ — + ^— = e 2 + ~V 2 2 + ^ + --^ + , 44 

2 pi 2 pi 2 ei 2 p 2 2 p 2 2 e 2 

where -Di( 2 ) is the electric induction in phase 1 (resp. 2), and €\ (resp. e 2 ) is the 
dielectric constant in phase 1 (resp. 2). We have neglected the magnetic contribution 
as the phase transition is mainly associated with a transfer of electric charges and is 
thus an electric effect. 

The appendix 2 has calculated the coupled elastic-electric modes in a piezoelectric 
medium. As an order of magnitude, we can estimate the amplitude of the electric 
induction discontinuity at the shock from the piezoelectric equation (|28|) of the ap- 
pendix 1 with E = (negligible external electric field) : D = da. Note that the 
discontinuity of the electric induction reflects the abrupt change of the polarization 
between the two phases. This gives the following equation which replaces (|23|) : 

e(p 2 ) + ^(4-3^ + 2) + -G.(^-D= C (p 1 ) + ^(4-3^+2) + ^G 2 (^-l) . 
2p 2 pi P2 ei Pi 2pi pf pi e 2 p 2 

(45) 

The discontinuity in stress and strain at the shock phase transition produces 
an impulse both mechanical and electrical. This generates seismic radiations and 
electric signals. This provides a natural mechanism for the generation of electric 
signals associated to earthquake rupture. The electric signal leads to a polarization 
wave which propagates in the crust and can be detected on the surface in favorable 
conditions. 

One should distinguish between two types of electric signals: 1) electromagnetic 
signals of very low frequency propagate as electromagnetic waves ; 2) polarization 
currents are carried along conducting paths within the crust. The first phenomenon 
is a fast process that occurs essentially with the explosive transformation while the 
second one is a slow phenomenon and can be detected after the conduction has 
brought the impulsive charge to the detectors. 

We now calculate the amplitude of the first electromagnetic process. For this, 
we use the wave equation (^,^) of the appendix 2, writing it in 3D, neglecting the 
small term ^d 2 << 1 and replacing the coupling term dGp-^^ by dp^- which acts 
as a source term corresponding to the passage of the shock, with stress discontinuity 
<7 2 — <J\- We thus have % = <J2 Ai " 1 , where At is the duration of the shock, i.e. the 
time it takes for the shock to pass over a point. In the inviscid limit, the width of 
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the shock shrinks to zero and At — > 0, thus leading to |j| — > {a 2 — cri)<5(t — t*(r)), 
where t* (r) is the time of arrival of the shock at f. We thus get finally the following 
wave equation with a source term due to the shock : 

&E-\^ = d^a 2 -a 1 )5'(t-t*) , (46) 

where A is the Laplacian and 5' is the derivative of the Dirac function. The formal 
solution of this equation is obtained by the method of Green function and reads : 



5' 



|r-f| [f*| 



ad r V c TJ , 

E(f, t) = p (a 2 - a,) / V ^ ^ . (47) 

47r J5 r — r 1 



We have replaced £* by K^-, which is the time of arrival of the shock at position f 
on the fault. The integral is carried over the surface covered by the shock. Since 
U < c (U is comparable to a sound wave velocity, while c is the speed of light in the 
medium), we can carry out the integration with the dirac function expressed in the 
r 1 variable and obtain in ID : 

This gives a signal with a long tail whose peak propagates at the shock velocity. 
The duration of the electric signal is thus directly proportional to the length of the 
rupture, and thus gives a direct information of the size of the event. 

In 3D, we define Ox as the direction along the long axis of the rupture, Oy is the 
direction parallel to the rupture and perpendicular to Ox and Oz is perperdicular to 
the rupture plane. Solving the integral in (|47|) with the dirac function, we get 



E(f,t) = ^(a 2 -a 1 ) jdy' % —^l , (49) 
47r J f — r \ 



where x'(y) is solution of 

U 
c 



Ut-x' = - y/(x - x') 2 + {y- y') 2 + z 2 , (50) 



and r' is expressed at the point of coordinates (x'(y'), y', 0). We notice that 
dy'^^ = dx'%^f^. We thus obtain 

a \f-~r 2 dx' ax' 



E{r,t) = — <7 2 - o-i) — \ B — — =- - — (a 2 - a l ) / dx dx , 51 
47r ax' \r — r'b\ 47r J r — r" 



where the index B refers to the contribution of the boundaries, denotes the position 
of the beginning and end of the rupture. There is thus a specific component of the 
signal radiated from the two edges of the rupture. 
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Figure captions : 



Figure 1: Starting from a double well configuration with Q_ (undeformed stable 
mineral) more stable than Q + (metastable mineral), the deformation applied to the 
Q_ phase creates a higher energy state which eventually becomes comparable to Q + 
(this is the effect of inclusion of dislocations, for instance). As a consequence, the 
system transforms into the metastable Q + . As the strain continues to increase, the 
free energy landscape deforms until a point where Q + becomes unstable and the 
mineral transform back into undeformed Q- . 

Figure 2: A one-dimensional chain is made of blocks linked to each other by 
energetic links which, when stressed beyond a given deformation threshold, rupture 
by releasing a burst of energy converted into kinetic energy transmitted to the blocks. 
The figure shows two successive bond ruptures that lead to velocity boosts to the 
ejected fragments on the left and to the boundary blocks. 
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